Circulatory levels of lysophosphatidylcholine species in obese adolescents: Findings from cross-sectional and prospective lipidomics analyses.

BACKGROUND AND AIMS: Obesity has reached epidemic proportions, emphasizing the importance of reliable biomarkers for detecting early metabolic alterations and enabling early preventative interventions. However, our understanding of the molecular mechanisms and specific lipid species associated with childhood obesity remains limited. Therefore, the aim of this study was to investigate plasma lipidomic signatures as potential biomarkers for adolescent obesity. METHODS AND RESULTS: A total of 103 individuals comprising overweight/obese (n = 46) and normal weight (n = 57) were randomly chosen from the baseline ORANGE (Obesity Reduction and Noncommunicable Disease Awareness through Group Education) cohort, having been followed up for a median of 7.1 years. Plasma lipidomic profiling was performed using the UHPLC-HRMS method. We used three different models adjusted for clinical covariates to analyze the data. Clustering methods were used to define metabotypes, which allowed for the stratification of subjects into subgroups with similar clinical and metabolic profiles. We observed that lysophosphatidylcholine (LPC) species like LPC.16.0, LPC.18.3, LPC.18.1, and LPC.20.3 were significantly (p < 0.05) associated with baseline and follow-up BMI in adolescent obesity. The association of LPC species with BMI remained consistently significant even after adjusting for potential confounders. Moreover, applying metabotyping using hierarchical clustering provided insights into the metabolic heterogeneity within the normal and obese groups, distinguishing metabolically healthy individuals from those with unhealthy metabolic profiles. CONCLUSION: The specific LPC levels were found to be altered and increased in childhood obesity, particularly during the follow-up. These findings suggest that LPC species hold promise as potential biomarkers of obesity in adolescents, including healthy and unhealthy metabolic profiles.

Genetic dissection of serum vaspin highlights its causal role in lipid metabolism.

OBJECTIVE: Vaspin (visceral adipose tissue derived serine protease inhibitor, SERPINA12) is associated with obesity-related metabolic traits, but its causative role is still elusive. The role of genetics in serum vaspin variability to establish its causal relationship with metabolically relevant traits was investigated. METHODS: A meta-analysis of genome-wide association studies for serum vaspin from six independent cohorts (N = 7446) was conducted. Potential functional variants of vaspin were included in Mendelian randomization (MR) analyses to assess possible causal pathways between vaspin and homeostasis model assessment and lipid traits. To further validate the MR analyses, data from Genotype-Tissue Expression (GTEx) were analyzed, db/db mice were treated with vaspin, and serum lipids were measured. RESULTS: A total of 468 genetic variants represented by five independent variants (rs7141073, rs1956709, rs4905216, rs61978267, rs73338689) within the vaspin locus were associated with serum vaspin (all p < 5×10-8 , explained variance 16.8%). MR analyses revealed causal relationships between serum vaspin and triglycerides, low-density lipoprotein, and total cholesterol. Gene expression correlation analyses suggested that genes, highly correlated with vaspin expression in adipose tissue, are enriched in lipid metabolic processes. Finally, in vivo vaspin treatment reduced serum triglycerides in obese db/db mice. CONCLUSIONS: The data show that serum vaspin is strongly determined by genetic variants within vaspin, which further highlight vaspin's causal role in lipid metabolism.

A Fecal Metabolite Signature of Impaired Fasting Glucose: Results From Two Independent Population-Based Cohorts.

Prediabetes is a metabolic condition associated with gut microbiome composition, although mechanisms remain elusive. We searched for fecal metabolites, a readout of gut microbiome function, associated with impaired fasting glucose (IFG) in 142 individuals with IFG and 1,105 healthy individuals from the UK Adult Twin Registry (TwinsUK). We used the Cooperative Health Research in the Region of Augsburg (KORA) cohort (318 IFG individuals, 689 healthy individuals) to replicate our findings. We linearly combined eight IFG-positively associated metabolites (1-methylxantine, nicotinate, glucuronate, uridine, cholesterol, serine, caffeine, and protoporphyrin IX) into an IFG-metabolite score, which was significantly associated with higher odds ratios (ORs) for IFG (TwinsUK: OR 3.9 [95% CI 3.02-5.02], P < 0.0001, KORA: OR 1.3 [95% CI 1.16-1.52], P < 0.0001) and incident type 2 diabetes (T2D; TwinsUK: hazard ratio 4 [95% CI 1.97-8], P = 0.0002). Although these are host-produced metabolites, we found that the gut microbiome is strongly associated with their fecal levels (area under the curve >70%). Abundances of Faecalibacillus intestinalis, Dorea formicigenerans, Ruminococcus torques, and Dorea sp. AF24-7LB were positively associated with IFG, and such associations were partially mediated by 1-methylxanthine and nicotinate (variance accounted for mean 14.4% [SD 5.1], P < 0.05). Our results suggest that the gut microbiome is linked to prediabetes not only via the production of microbial metabolites but also by affecting intestinal absorption/excretion of host-produced metabolites and xenobiotics, which are correlated with the risk of IFG. Fecal metabolites enable modeling of another mechanism of gut microbiome effect on prediabetes and T2D onset. ARTICLE HIGHLIGHTS: Prediabetes is a metabolic condition associated with gut microbiome composition, although mechanisms remain elusive. We investigated whether there is a fecal metabolite signature of impaired fasting glucose (IFG) and the possible underlying mechanisms of action. We identified a fecal metabolite signature of IFG associated with prevalent IFG in two independent cohorts and incident type 2 diabetes in a subanalysis. Although the signature consists of metabolites of nonmicrobial origin, it is strongly correlated with gut microbiome composition. Fecal metabolites enable modeling of another mechanism of gut microbiome effect on prediabetes by affecting intestinal absorption or excretion of host compounds and xenobiotics.

GWAS of random glucose in 476,326 individuals provide insights into diabetes pathophysiology, complications and treatment stratification.

Conventional measurements of fasting and postprandial blood glucose levels investigated in genome-wide association studies (GWAS) cannot capture the effects of DNA variability on 'around the clock' glucoregulatory processes. Here we show that GWAS meta-analysis of glucose measurements under nonstandardized conditions (random glucose (RG)) in 476,326 individuals of diverse ancestries and without diabetes enables locus discovery and innovative pathophysiological observations. We discovered 120 RG loci represented by 150 distinct signals, including 13 with sex-dimorphic effects, two cross-ancestry and seven rare frequency signals. Of these, 44 loci are new for glycemic traits. Regulatory, glycosylation and metagenomic annotations highlight ileum and colon tissues, indicating an underappreciated role of the gastrointestinal tract in controlling blood glucose. Functional follow-up and molecular dynamics simulations of lower frequency coding variants in glucagon-like peptide-1 receptor (GLP1R), a type 2 diabetes treatment target, reveal that optimal selection of GLP-1R agonist therapy will benefit from tailored genetic stratification. We also provide evidence from Mendelian randomization that lung function is modulated by blood glucose and that pulmonary dysfunction is a diabetes complication. Our investigation yields new insights into the biology of glucose regulation, diabetes complications and pathways for treatment stratification.

Genetics of circulating inflammatory proteins identifies drivers of immune-mediated disease risk and therapeutic targets.

Circulating proteins have important functions in inflammation and a broad range of diseases. To identify genetic influences on inflammation-related proteins, we conducted a genome-wide protein quantitative trait locus (pQTL) study of 91 plasma proteins measured using the Olink Target platform in 14,824 participants. We identified 180 pQTLs (59 cis, 121 trans). Integration of pQTL data with eQTL and disease genome-wide association studies provided insight into pathogenesis, implicating lymphotoxin-α in multiple sclerosis. Using Mendelian randomization (MR) to assess causality in disease etiology, we identified both shared and distinct effects of specific proteins across immune-mediated diseases, including directionally discordant effects of CD40 on risk of rheumatoid arthritis versus multiple sclerosis and inflammatory bowel disease. MR implicated CXCL5 in the etiology of ulcerative colitis (UC) and we show elevated gut CXCL5 transcript expression in patients with UC. These results identify targets of existing drugs and provide a powerful resource to facilitate future drug target prioritization.

Genetic analysis of blood molecular phenotypes reveals common properties in the regulatory networks affecting complex traits.

We evaluate the shared genetic regulation of mRNA molecules, proteins and metabolites derived from whole blood from 3029 human donors. We find abundant allelic heterogeneity, where multiple variants regulate a particular molecular phenotype, and pleiotropy, where a single variant associates with multiple molecular phenotypes over multiple genomic regions. The highest proportion of share genetic regulation is detected between gene expression and proteins (66.6%), with a further median shared genetic associations across 49 different tissues of 78.3% and 62.4% between plasma proteins and gene expression. We represent the genetic and molecular associations in networks including 2828 known GWAS variants, showing that GWAS variants are more often connected to gene expression in trans than other molecular phenotypes in the network. Our work provides a roadmap to understanding molecular networks and deriving the underlying mechanism of action of GWAS variants using different molecular phenotypes in an accessible tissue.

A polygenic and family risk score are both independently associated with risk of type 2 diabetes in a population-based study.

The availability of polygenic scores for type 2 diabetes (T2D) raises the question, whether assessing family history might become redundant. However, family history not only involves shared genetics, but also shared environment. It was the aim of this study to assess the independent and combined effects of one family risk score (FamRS) and a polygenic score (PGS) on prevalent and incident T2D risk in a population-based study from Germany (n = 3071). The study was conducted in 2004/2005 with up to 12 years of follow-up. The FamRS takes into account not only the number of diseased first grade relatives, but also age at onset of the relatives and age of participants. 256 prevalent and additional 163 incident T2D cases were recorded. Prevalence of T2D increased sharply for those within the top quantile of the PGS distribution resulting in an OR of 19.16 (p < 2 × 10-16) for the top 20% compared to the remainder of the population, independent of age, sex, BMI, physical activity and FamRS. On the other hand, having a very strong family risk compared to average was still associated with an OR of 2.78 (p = 0.001), independent of the aforementioned factors and the PGS. The PGS and FamRS were only slightly correlated (r2Spearman = 0.018). The combined contribution of both factors varied with varying age-groups, though, with decreasing influence of the PGS with increasing age. To conclude, both, genetic information and family history are relevant for the prediction of T2D risk and might be used for identification of high risk groups to personalize prevention measures.

Metabolic Signatures Elucidate the Effect of Body Mass Index on Type 2 Diabetes.

Obesity plays an important role in the development of insulin resistance and diabetes, but the molecular mechanism that links obesity and diabetes is still not completely understood. Here, we used 146 targeted metabolomic profiles from the German KORA FF4 cohort consisting of 1715 participants and associated them with obesity and type 2 diabetes. In the basic model, 83 and 51 metabolites were significantly associated with body mass index (BMI) and T2D, respectively. Those metabolites are branched-chain amino acids, acylcarnitines, lysophospholipids, or phosphatidylcholines. In the full model, 42 and 3 metabolites were significantly associated with BMI and T2D, respectively, and replicate findings in the previous studies. Sobel mediation testing suggests that the effect of BMI on T2D might be mediated via lipids such as sphingomyelin (SM) C16:1, SM C18:1 and diacylphosphatidylcholine (PC aa) C38:3. Moreover, mendelian randomization suggests a causal relationship that BMI causes the change of SM C16:1 and PC aa C38:3, and the change of SM C16:1, SM C18:1, and PC aa C38:3 contribute to T2D incident. Biological pathway analysis in combination with genetics and mice experiments indicate that downregulation of sphingolipid or upregulation of phosphatidylcholine metabolism is a causal factor in early-stage T2D pathophysiology. Our findings indicate that metabolites like SM C16:1, SM C18:1, and PC aa C38:3 mediate the effect of BMI on T2D and elucidate their role in obesity related T2D pathologies.

Loci for insulin processing and secretion provide insight into type 2 diabetes risk.

Insulin secretion is critical for glucose homeostasis, and increased levels of the precursor proinsulin relative to insulin indicate pancreatic islet beta-cell stress and insufficient insulin secretory capacity in the setting of insulin resistance. We conducted meta-analyses of genome-wide association results for fasting proinsulin from 16 European-ancestry studies in 45,861 individuals. We found 36 independent signals at 30 loci (p value < 5 × 10-8), which validated 12 previously reported loci for proinsulin and ten additional loci previously identified for another glycemic trait. Half of the alleles associated with higher proinsulin showed higher rather than lower effects on glucose levels, corresponding to different mechanisms. Proinsulin loci included genes that affect prohormone convertases, beta-cell dysfunction, vesicle trafficking, beta-cell transcriptional regulation, and lysosomes/autophagy processes. We colocalized 11 proinsulin signals with islet expression quantitative trait locus (eQTL) data, suggesting candidate genes, including ARSG, WIPI1, SLC7A14, and SIX3. The NKX6-3/ANK1 proinsulin signal colocalized with a T2D signal and an adipose ANK1 eQTL signal but not the islet NKX6-3 eQTL. Signals were enriched for islet enhancers, and we showed a plausible islet regulatory mechanism for the lead signal in the MADD locus. These results show how detailed genetic studies of an intermediate phenotype can elucidate mechanisms that may predispose one to disease.

Epigenome-Wide Association Study Reveals CpG Sites Associated with Thyroid Function and Regulatory Effects on KLF9.

Background: Thyroid hormones play a key role in differentiation and metabolism and are known regulators of gene expression through both genomic and epigenetic processes including DNA methylation. The aim of this study was to examine associations between thyroid hormones and DNA methylation. Methods: We carried out a fixed-effect meta-analysis of epigenome-wide association study (EWAS) of blood DNA methylation sites from 8 cohorts from the ThyroidOmics Consortium, incorporating up to 7073 participants of both European and African ancestry, implementing a discovery and replication stage. Statistical analyses were conducted using normalized beta CpG values as dependent and log-transformed thyrotropin (TSH), free thyroxine, and free triiodothyronine levels, respectively, as independent variable in a linear model. The replicated findings were correlated with gene expression levels in whole blood and tested for causal influence of TSH and free thyroxine by two-sample Mendelian randomization (MR). Results: Epigenome-wide significant associations (p-value <1.1E-7) of three CpGs for free thyroxine, five for free triiodothyronine, and two for TSH concentrations were discovered and replicated (combined p-values = 1.5E-9 to 4.3E-28). The associations included CpG sites annotated to KLF9 (cg00049440) and DOT1L (cg04173586) that overlap with all three traits, consistent with hypothalamic-pituitary-thyroid axis physiology. Significant associations were also found for CpGs in FKBP5 for free thyroxine, and at CSNK1D/LINCO1970 and LRRC8D for free triiodothyronine. MR analyses supported a causal effect of thyroid status on DNA methylation of KLF9. DNA methylation of cg00049440 in KLF9 was inversely correlated with KLF9 gene expression in blood. The CpG at CSNK1D/LINC01970 overlapped with thyroid hormone receptor alpha binding peaks in liver cells. The total additive heritability of the methylation levels of the six significant CpG sites was between 25% and 57%. Significant methylation QTLs were identified for CpGs at KLF9, FKBP5, LRRC8D, and CSNK1D/LINC01970. Conclusions: We report novel associations between TSH, thyroid hormones, and blood-based DNA methylation. This study advances our understanding of thyroid hormone action particularly related to KLF9 and serves as a proof-of-concept that integrations of EWAS with other -omics data can provide a valuable tool for unraveling thyroid hormone signaling in humans by complementing and feeding classical in vitro and animal studies.

Differences and commonalities in the genetic architecture of protein quantitative trait loci in European and Arab populations.

Polygenic scores (PGS) can identify individuals at risk of adverse health events and guide genetics-based personalized medicine. However, it is not clear how well PGS translate between different populations, limiting their application to well-studied ethnicities. Proteins are intermediate traits linking genetic predisposition and environmental factors to disease, with numerous blood circulating protein levels representing functional readouts of disease-related processes. We hypothesized that studying the genetic architecture of a comprehensive set of blood-circulating proteins between a European and an Arab population could shed fresh light on the translatability of PGS to understudied populations. We therefore conducted a genome-wide association study with whole-genome sequencing data using 1301 proteins measured on the SOMAscan aptamer-based affinity proteomics platform in 2935 samples of Qatar Biobank and evaluated the replication of protein quantitative traits (pQTLs) from European studies in an Arab population. Then, we investigated the colocalization of shared pQTL signals between the two populations. Finally, we compared the performance of protein PGS derived from a Caucasian population in a European and an Arab cohort. We found that the majority of shared pQTL signals (81.8%) colocalized between both populations. About one-third of the genetic protein heritability was explained by protein PGS derived from a European cohort, with protein PGS performing ~20% better in Europeans when compared to Arabs. Our results are relevant for the translation of PGS to non-Caucasian populations, as well as for future efforts to extend genetic research to understudied populations.

Network reconstruction for trans acting genetic loci using multi-omics data and prior information.

BACKGROUND: Molecular measurements of the genome, the transcriptome, and the epigenome, often termed multi-omics data, provide an in-depth view on biological systems and their integration is crucial for gaining insights in complex regulatory processes. These data can be used to explain disease related genetic variants by linking them to intermediate molecular traits (quantitative trait loci, QTL). Molecular networks regulating cellular processes leave footprints in QTL results as so-called trans-QTL hotspots. Reconstructing these networks is a complex endeavor and use of biological prior information can improve network inference. However, previous efforts were limited in the types of priors used or have only been applied to model systems. In this study, we reconstruct the regulatory networks underlying trans-QTL hotspots using human cohort data and data-driven prior information. METHODS: We devised a new strategy to integrate QTL with human population scale multi-omics data. State-of-the art network inference methods including BDgraph and glasso were applied to these data. Comprehensive prior information to guide network inference was manually curated from large-scale biological databases. The inference approach was extensively benchmarked using simulated data and cross-cohort replication analyses. Best performing methods were subsequently applied to real-world human cohort data. RESULTS: Our benchmarks showed that prior-based strategies outperform methods without prior information in simulated data and show better replication across datasets. Application of our approach to human cohort data highlighted two novel regulatory networks related to schizophrenia and lean body mass for which we generated novel functional hypotheses. CONCLUSIONS: We demonstrate that existing biological knowledge can improve the integrative analysis of networks underlying trans associations and generate novel hypotheses about regulatory mechanisms.

Fecal Bile Acids and Neutral Sterols Are Associated with Latent Microbial Subgroups in the Human Gut.

Bile acids, neutral sterols, and the gut microbiome are intricately intertwined and each affects human health and metabolism. However, much is still unknown about this relationship. This analysis included 1280 participants of the KORA FF4 study. Fecal metabolites (primary and secondary bile acids, plant and animal sterols) were analyzed using a metabolomics approach. Dirichlet regression models were used to evaluate associations between the metabolites and twenty microbial subgroups that were previously identified using latent Dirichlet allocation. Significant associations were identified between 12 of 17 primary and secondary bile acids and several of the microbial subgroups. Three subgroups showed largely positive significant associations with bile acids, and six subgroups showed mostly inverse associations with fecal bile acids. We identified a trend where microbial subgroups that were previously associated with "healthy" factors were here inversely associated with fecal bile acid levels. Conversely, subgroups that were previously associated with "unhealthy" factors were positively associated with fecal bile acid levels. These results indicate that further research is necessary regarding bile acids and microbiota composition, particularly in relation to metabolic health.

A randomization-based causal inference framework for uncovering environmental exposure effects on human gut microbiota.

Statistical analysis of microbial genomic data within epidemiological cohort studies holds the promise to assess the influence of environmental exposures on both the host and the host-associated microbiome. However, the observational character of prospective cohort data and the intricate characteristics of microbiome data make it challenging to discover causal associations between environment and microbiome. Here, we introduce a causal inference framework based on the Rubin Causal Model that can help scientists to investigate such environment-host microbiome relationships, to capitalize on existing, possibly powerful, test statistics, and test plausible sharp null hypotheses. Using data from the German KORA cohort study, we illustrate our framework by designing two hypothetical randomized experiments with interventions of (i) air pollution reduction and (ii) smoking prevention. We study the effects of these interventions on the human gut microbiome by testing shifts in microbial diversity, changes in individual microbial abundances, and microbial network wiring between groups of matched subjects via randomization-based inference. In the smoking prevention scenario, we identify a small interconnected group of taxa worth further scrutiny, including Christensenellaceae and Ruminococcaceae genera, that have been previously associated with blood metabolite changes. These findings demonstrate that our framework may uncover potentially causal links between environmental exposure and the gut microbiome from observational data. We anticipate the present statistical framework to be a good starting point for further discoveries on the role of the gut microbiome in environmental health.

Association between a polygenic and family risk score on the prevalence and incidence of myocardial infarction in the KORA-F3 study.

BACKGROUND AND AIMS: Genetic risk scores for common diseases as myocardial infarction (MI) gain increasing attention for individual's risk prediction. One might wonder if assessing family history becomes redundant. It was the aim of this study to evaluate the amount of shared information between family history and genetic risk scores and to assess their independent and combined effects on prevalent and incident MI risk. METHODS: A genome wide polygenic risk score (PGS) and one family risk score (FamRS) were calculated in a population-based study from Southern Germany (n = 3071) with up to 11 years of follow-up. Logistic and Cox Regression models were used adjusting for lifestyle and classical risk factors. RESULTS: A right shift in MI risk for increasing values of PGS was found, with. considerably increasing ORs along the top quantiles of PGS (OR = 3.03 for top 10%; OR = 5.55 for top 2.5%). The PGS was not associated with incident MI cases, though. The FamRS was significantly associated with both prevalent and incident MI cases with an OR of 2.9 for participants with a strong positive family history compared to average. ORs and HRs did hardly change in a combined model including both measures, indicating independent contribution to MI risk. The simultaneous addition of PGS and FamRS to a model including classical risk factors significantly enhanced prediction for prevalent cases and non-cases (p = 3.28 × 10-5). CONCLUSIONS: These findings emphasize that both genetic information and family history are relevant for the determination of MI risk and that neither of them can replace the other.

A Family and a Genome-Wide Polygenic Risk Score Are Independently Associated With Stroke in a Population-Based Study.

BACKGROUND: Positive family history and genetic risk scores have been shown to independently capture those individuals with high risk for stroke. The aim of our study was to evaluate the amount of shared information between family history and genetic risk and to investigate their combined effect on the association with prevalent and incident stroke cases. METHODS: We obtained a family risk score (FamRS), weighted for disease onset and family size as well as genome-wide polygenic risk score (PGS) including over 3.2 million single-nucleotide polymorphisms in the population-based prospective KORA F3 (Cooperative Health Research in the Region of Augsburg) study (n=3071) from Southern Germany. FamRS and PGS were evaluated separately and combined. The measures were once treated as continuous variables but also divided in the highest 20%, 10%, 5%, and 1% percentiles. Odds ratios via logistic regression and hazard ratios via Cox regression were estimated. A stroke event was defined as a hospitalization for stroke that was self-reported in a standardized interview by certified and supervised personnel. RESULTS: The FamRS outperformed other simplified family measures such as affected parents or number of affected family members. FamRS and PGS were not correlated, and no individuals were observed with both very high FamRS and very high PGS (top 1% percentile). In a combined model, both FamRS and PGS were independently from each other associated with risk of stroke, also independent of other traditional risk factors (p [FamRS]=0.02, p [PGS]=0.005). Individuals in the top 1% of either FamRS or PGS were found to have >5-fold risk for stroke (odds ratios, 5.82 [95% CI, 2.08-14]; P=0.0002). The results for incident stroke events showed the same trend but were not significant. CONCLUSIONS: Our study shows that a family risk score and PGS capture different information concerning individual stroke risk. Combining the risk measures FamRS and PGS increases predictive power, as demonstrated in a population-based study.

Colocalization analysis of pancreas eQTLs with risk loci from alcoholic and novel non-alcoholic chronic pancreatitis GWAS suggests potential disease causing mechanisms.

BACKGROUND: Previous genome-wide association studies (GWAS) identified genome-wide significant risk loci in chronic pancreatitis and investigated underlying disease causing mechanisms by simple overlaps with expression quantitative trait loci (eQTLs), a procedure which may often result in false positive conclusions. METHODS: We conducted a GWAS in 584 non-alcoholic chronic pancreatitis (NACP) patients and 6040 healthy controls. Next, we applied Bayesian colocalization analysis of identified genome-wide significant risk loci from both, our recently published alcoholic chronic pancreatitis (ACP) and the novel NACP dataset, with pancreas eQTLs from the GTEx V8 European cohort to prioritize candidate causal genes and extracted credible sets of shared causal variants. RESULTS: Variants at the CTRC (p = 1.22 × 10-21) and SPINK1 (p = 6.59 × 10-47) risk loci reached genome-wide significance in NACP. CTRC risk variants colocalized with CTRC eQTLs in ACP (PP4 = 0.99, PP4/PP3 = 95.51) and NACP (PP4 = 0.99, PP4/PP3 = 95.46). For both diseases, the 95% credible set of shared causal variants consisted of rs497078 and rs545634. CLDN2-MORC4 risk variants colocalized with CLDN2 eQTLs in ACP (PP4 = 0.98, PP4/PP3 = 42.20) and NACP (PP4 = 0.67, PP4/PP3 = 7.18), probably driven by the shared causal variant rs12688220. CONCLUSIONS: A shared causal CTRC risk variant might unfold its pathogenic effect in ACP and NACP by reducing CTRC expression, while the CLDN2-MORC4 shared causal variant rs12688220 may modify ACP and NACP risk by increasing CLDN2 expression.

Epigenome-wide association study of incident type 2 diabetes: a meta-analysis of five prospective European cohorts.

AIMS/HYPOTHESIS: Type 2 diabetes is a complex metabolic disease with increasing prevalence worldwide. Improving the prediction of incident type 2 diabetes using epigenetic markers could help tailor prevention efforts to those at the highest risk. The aim of this study was to identify predictive methylation markers for incident type 2 diabetes by combining epigenome-wide association study (EWAS) results from five prospective European cohorts. METHODS: We conducted a meta-analysis of EWASs in blood collected 7-10 years prior to type 2 diabetes diagnosis. DNA methylation was measured with Illumina Infinium Methylation arrays. A total of 1250 cases and 1950 controls from five longitudinal cohorts were included: Doetinchem, ESTHER, KORA1, KORA2 and EPIC-Norfolk. Associations between DNA methylation and incident type 2 diabetes were examined using robust linear regression with adjustment for potential confounders. Inverse-variance fixed-effects meta-analysis of cohort-level individual CpG EWAS estimates was performed using METAL. The methylGSA R package was used for gene set enrichment analysis. Confirmation of genome-wide significant CpG sites was performed in a cohort of Indian Asians (LOLIPOP, UK). RESULTS: The meta-analysis identified 76 CpG sites that were differentially methylated in individuals with incident type 2 diabetes compared with control individuals (p values <1.1 × 10-7). Sixty-four out of 76 (84.2%) CpG sites were confirmed by directionally consistent effects and p values <0.05 in an independent cohort of Indian Asians. However, on adjustment for baseline BMI only four CpG sites remained genome-wide significant, and addition of the 76 CpG methylation risk score to a prediction model including established predictors of type 2 diabetes (age, sex, BMI and HbA1c) showed no improvement (AUC 0.757 vs 0.753). Gene set enrichment analysis of the full epigenome-wide results clearly showed enrichment of processes linked to insulin signalling, lipid homeostasis and inflammation. CONCLUSIONS/INTERPRETATION: By combining results from five European cohorts, and thus significantly increasing study sample size, we identified 76 CpG sites associated with incident type 2 diabetes. Replication of 64 CpGs in an independent cohort of Indian Asians suggests that the association between DNA methylation levels and incident type 2 diabetes is robust and independent of ethnicity. Our data also indicate that BMI partly explains the association between DNA methylation and incident type 2 diabetes. Further studies are required to elucidate the underlying biological mechanisms and to determine potential causal roles of the differentially methylated CpG sites in type 2 diabetes development.

Genetic variation influencing DNA methylation provides insights into molecular mechanisms regulating genomic function.

We determined the relationships between DNA sequence variation and DNA methylation using blood samples from 3,799 Europeans and 3,195 South Asians. We identify 11,165,559 SNP-CpG associations (methylation quantitative trait loci (meQTL), P < 10-14), including 467,915 meQTL that operate in trans. The meQTL are enriched for functionally relevant characteristics, including shared chromatin state, High-throuhgput chromosome conformation interaction, and association with gene expression, metabolic variation and clinical traits. We use molecular interaction and colocalization analyses to identify multiple nuclear regulatory pathways linking meQTL loci to phenotypic variation, including UBASH3B (body mass index), NFKBIE (rheumatoid arthritis), MGA (blood pressure) and COMMD7 (white cell counts). For rs6511961 , chromatin immunoprecipitation followed by sequencing (ChIP-seq) validates zinc finger protein (ZNF)333 as the likely trans acting effector protein. Finally, we used interaction analyses to identify population- and lineage-specific meQTL, including rs174548 in FADS1, with the strongest effect in CD8+ T cells, thus linking fatty acid metabolism with immune dysregulation and asthma. Our study advances understanding of the potential pathways linking genetic variation to human phenotype.